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Abstract. A signature independent formalism is created and utilized to determine the ge- 
neral second-order symmetry operators for Dirac's equation on two-dimensional Lorentzian 
I spin manifolds. The formalism is used to characterize the ortho normal frames and metrics 

that permit the solution of Dirac's equation by separation of variables in the case where 
a second-order symmetry operator underlies the separation. Separation of variables in com- 
■ plex variables on two-dimensional Minkowski space is also considered. 
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1 Introduction 



■ This paper is a contribution to the study of the separability theory for Dirac's equation to which 

Professor Miller has made important contributions [171 [TSj [24] . Exact solutions to Dirac's re- 
lativistic wave equation by means of the method of separation of variables have been studied 
since the equation was postulated in 1928. Indeed, the solution for the hydrogen atom may be 
obtained by this method. While there is a well developed theory of separation of variables for the 
Hamilton-Jacobi equation, and the Schrodinger equation based on the existence of valence two 
characteristic Killing tensors which define respectively quadratic first integrals and second-order 
symmetry operators for these equations (see [231 [T6t [H [12]) an analogous theory for the Dirac 
equation is still in its early stages. The complications arise from the fact that one is dealing 
with a system of first-order partial differential equations whose derivation from the invariant 
Dirac equation depends not only on the choice of coordinate system but also on the choice of an 
orthonormal moving frame and representation for the Dirac matrices with respect to which the 
components of the unknown spinor are defined. Further complications arise if the background 
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space-time is assumed to be curved. Much of the progress in the theory has been stimulated by 
developments in Einstein's general theory of relativity where one studies first quantized relativis- 
tic electrons on curved background space-times of physical interest such as the Schwarzschild 
and Kerr black hole space-times. This work required the preliminary development of a theory 
of spinors on general pseudo-Riemmanian manifolds (see [9l \5\ [TOl [11] )• The solution of the 
Dirac equation in the Reissner-Nordstrom solution was apparently first obtained by Brill and 
Wheeler in 1957 [2] who separated the equations for the spinor components in standard orthogo- 
nal Schwarzschild coordinates with respect to a moving frame adapted to the coordinate curves. 
A comparable separable solution in the Kerr solution was found by Chandrasekhar in 1976 [6] by 
use of an ingenious separation ansatz involving Boyer-Lindquist coordinates and the Kinnersley 
tetrad. The separability property was characterized invariantly by Carter and McLenaghan 
in terms of a first-order differential operator constructed from the valence two Yano-Killing 
tensor that exists in the Kerr solution, that commutes with the Dirac operator and that admits 
the separable solutions as eigenfunctions with the separation constant as eigenvalue. Study of 
this example led Miller [24] to propose the theory of a factorizable system for first-order sys- 
tems of Dirac type in the context of which the separability property may be characterized by the 
existence of a certain system of commuting first-order symmetry operators. While this theory 
includes the Dirac equation on the Kerr solution and its generalizations [15] it is apparently not 
complete since as is shown by Fels and Kamran [15] there exist systems of the Dirac type whose 
separability is characterized by second-order symmetry operators. The work begun by these 
authors has been continued by Smith [25], Fatibene, McLenaghan and Smith [13], McLenaghan 
and Rastelli [20], and Fatibene, McLenaghan, Rastelli and Smith [13] who studied the problem in 
the simplest possible setting namely on two-dimensional Riemmanian spin manifolds. The mo- 
tivation for working in the lowest permitted dimension is that it is possible to examine in detail 
the different possible scenarios that arise from the separation ansatz and the imposition of the 
separation paradigm that the separation be characterized by a symmetry operator admitting the 
separable solutions as eigenfunctions. The insight obtained from this approach may help suggest 
an approach to take for the construction of a general separability theory for Dirac type equa- 
tions. Indeed in [20] systems of two first-order linear partial equations of Dirac type which admit 
multiplicative separation of variables in some arbitrary coordinate system and whose separation 
constants are associated with commuting differential operators are exhaustively characterized. 
The requirement that the original system arises from the Dirac equation on some two-dimensional 
Riemannian spin manifold allows the local characterization of the orthonormal frames and met- 
rics admitting separation of variables for the equation and the determination of the symmetry 
operators associated to the separation. The paper [13] takes this research in a different but 
closely related direction. Following earlier work of McLenaghan and Spindel [22] and Kamran 
and McLenaghan [19] where the first-order symmetry operators of the Dirac equation where com- 
puted on four-dimensional Lorentzian spin manifolds and McLenaghan, Smith and Walker [21] 
where the second-order operators were determined in terms of a two-component spinor formal- 
ism, the most general second-order linear differential operator which commutes with the Dirac 
operator on a general two-dimensional Riemannian spin manifold is obtained. Further it is shown 
that the operator is characterized in terms of Killing vectors and valence two Killing tensors 
defined on the background manifold. The derivation is manifestly covariant: the calculations 
are done in a general orthonormal frame without the choice of a particular set of Dirac matrices. 

The purpose of the present paper is to extend the results of [20] and [13] described above 
to the case of two-dimensional Lorentzian spin manifolds. One of the main achievements of 
the paper is the creation of a formalism that permits the simultaneous treatment of both the 
Riemannian and Lorentzian cases. Further, we extend the results of [8], where Hamilton-Jacobi 
separability separability is studied in complex variables on two-dimensional Minkowski space, 
to the Dirac equation. 
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The paper is organized as follows. In Section [2] we summarize the basic properties of two- 
dimensional spin manifolds required for the subsequent calculations. Section [3] is devoted to 
the derivation of the form of the general second-order linear differential operator which com- 
mutes with the Dirac operator. We show that this operator is characterized by a valence two 
Killing tensor field, two Killing vector fields and two scalar fields defined on the background 
spin-manifold. In Section |4] we develop a formalism based on [20] which enables us to study si- 
multaneously separation of variables for the Dirac equation on both Riemannian and Lorentzian 
spin manifolds. All possible cases where separation occurs are determined. In Section [5] we es- 
tablish a link between the second-order symmetry operators obtained formally in Section [3] and 
the second-order symmetry operator underlying the separation of variables scheme considered in 
Section m In Section [6] separation of Dirac's equation in complex variables is considered on two- 
dimensional Minkowski space. Section [7] contains the appendix. The notation and conventions 
of this paper are consistent with |14j . 



2 Framework 

Let M be a connected, paracompact, two-dimensional spin manifold. Let us consider both the 
Euclidean r] = (2,0) and the Lorentzian r] = (1, 1) signatures. With an abuse of notation, let r] 
also denote the canonical form induced by the signature and the determinant of such quadratic 
form, namely one has r/ = ±1. We will keep this sign as an undetermined parameter in this 
paper, since we want to consider both cases at once. 

We know that a representation of the Clifford algebra is induced by a set of Dirac matrices 7a 
such that they satisfy the Dirac condition 

7a 7b + 767a = 2??afe^- 

The generators of the even Clifford algebra are /, -ja and 7 := 7172. Therefore the most general 
element of the Spin(r7) group is 

S = al + bj 

with + = 1. 

From the theory, we know that it is possible to define a covering map between Spin(r7) 
and S0(?7). Let be a generic element of SO(r/). Then 

Vab = llVcdt 
and the covering map in matrix form is 

^^^>~\-2i]ab a^-7]by- 

Let P — )• M be a suitable Spin(r7) principal bundle, such that it allows global maps A : P — )• 
L{M) of the spin bundle into the general frame bundle of M. The local expression of such maps 
is given by spin frames Ca- 

A spin frame induces a metric of signature 77 and the corresponding spin connection: 

where T^^ denotes the Levi-Civita connection of the induced metric g^y, which in turn induces 
the covariant derivative as = e^Va. 

It follows from our setting that the inner product r]ab will raise and lower Latin indices, while 
the Greek ones are raised and lowered by the induced metric g^j^u- 
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Finally, the covariant derivative of spinors is defined as 

v^^ := + r^V', (1) 

where 

The covariant derivative ([T|) is invariant under spin transformations in view of the following 
property 

A spin transformation is an automorphism of the spin bundle P. If we define the fiber co- 
ordinates of P as {x,S), the local expressions of the spin transformation are x'^ = x'^{x), 
S' = (j){x)S, where (j) : U ^ Spin(?7). 

As we know, it acts as a left group on spinors and spin frames (/^ = la{^)) 

where / : Spin(ry) — )• S0(t7) and Jy is the Jacobian of the coordinate transformation. Therefore 
the spin connection transforms as 

Looking back to the covariant derivative, one can prove that the commutator can be related to 
the curvature 

[V^,V.]V = hV'eafe<.. (2) 

The advantage of working in dimension two is that the Riemann tensor R^j^y has only one 
independent component that can be written as a function of the Ricci scalar R. Hence the 
identity ([2]) may be rewritten as 

[Vc,Vd]i^ = \-fi^e,dR. (3) 
Similarly the following property holds 

3 First- and second-order symmetry operators 

In our framework, the Dirac equation has the form 

Dt/; = i7"VaV' - mV- = 0. (5) 
An operator is a symmetry operator for the Dirac equation if 

[K, B] = 0. (6) 
The most general operator of the second-order has the form 

IK = ^f'^Vab + 1F"V^ + GI, 
where E'^^ F'^, G are algebraic matrix coefficients to be determined. 
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We aim to determine coefficients so that condition ^ holds true. Considering the Ricci's 
identities ^ and ([ID, we expand the symmetry equation ^ for the coefficients and we 
obtain 

E^^'^j"^ - -f^'^E"''^ = 0, f(''7*) - 7(*F'') = 7''VcE''^ 

7'VeG = f (F^^ + t'F'^) 7eac + (h"^"' + I^"V) h^bcVbR. (7) 

To write the system ([7]) we make use of a characterization of second- and third-order covariant 
derivatives which is shown in the Appendix. 

Let us begin with the first equation in ([7]). We know that the coefficients of E"'' are zero-order 
matrix operators (i.e. not differential) which can be expanded in the basis of C{r]) 

where the coefficients are point functions in M. Here we are using the fact that /, and 7 
form a basis, since they are linearly independent. Hence, the first equation can be rewritten as 

which can be solved to obtain 

where a'^ are the coefficients along the frame Ca := Cadf^ of an arbitrary vector field a. 
Applying these conditions, the operator E is of the form 

Now we consider the second condition of ([7]). As we did for E, we expand F in the basis 

F'^ = r/+/,v + r7, 

and substitute it back in the condition ([7|). We obtain 

V^a^^T]^^^ = 0, /i'^e^)" = V,a('^e^)^ -If'^^'e^l = V,e'^^ (8) 

The first equation says that a" is a Killing vector of g or possibly the zero vector. Through 
an explicit calculation, we prove that the most general tensor satisfying the second condition 
is A6^. Hence 

Finally, we substitute into the third equation of ([8]) obtaining 

-2f-e'le'., = V,e-'e'.,. (9) 
Taking the trace we get 

fb = IVee"^^, = iV^e-^^,,, (10) 

which shows that is uniquely determined by e""^. However Q contains six equations, of which 
only one has been used. The other five equations are exploited by back substituting pO|) into Q 
to obtain an equation for e'^^ alone, namely 
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This is an integrability condition for e"''' which may be written as 
V^'^e''^) = 0. 

This equation shows that e"* is a KilHng tensor of g. 
We shaU now consider the third condition of d?]) 

Considering the usual expansion G = gl + ga'j'^ + 97; we obtain the following system of 
equations 

= VaA - §0^1 + §0^1, 

^9b^.. — (,VfeA- -jabje.. — -jab^.. — -jo^b^- j 

O^A^f*^^,^ V7 fC' \ V c c byjeyjd^ac , R „ab , R^cd^a- ^ b 

The first equation means that ^4 E C, while the second implies gb = — ^a^. By splitting the 
third into its symmetric and antisymmetric parts, it results that 

g = hba^'r 

for the antisymmetric part, while the symmetric part can be expanded to obtain 
V('^(/) - V^e^y) = 0. 

It follows that there exists a Killing vector or a zero- vector such that Q"" = — Vcc"'^ from 
which we obtain 

We summarize all the results obtained so far 
gab ^ e«6j + 2a(V), 

F'^ = {C + Vee"^)/ + {^"Vca" + A^") + 1 e^^V^e^^T, 
where 

^ G C, V^'^e"'') = 0, V^'^a*) = 0, V^'^C^^ = 0. 
It remains to consider the fourth equation in ([7]) 

T^VeG = f (E^^ + t'F'^) leac + (h'^E"' + E^V) l^ebc^^/bR. 

which implies the only additional condition: 

= -\Vb{Re'''). (11) 

This equation locally determines g if and only if the right hand side is a closed 1-form. We 
call (|lip the integrability condition. 

We now focus on finding condition for first-order operators. 

We can easily obtain the conditions by setting to zero e"^ and a". In particular (llip is trivially 
satisfied and 5 G C. We thus obtain that the most general first-order symmetry operator may 
be written as 

^ab ^ jpa ^ ^aj ^ ^^a ^ G = gl + ie^.V^C^, 

where A,geC and V^'^C^^ = 0. 
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4 Separation of variables 

Let us start with the Dirac condition 

^^^b ^ ^b^a ^ 2if^I. (12) 

We fix the metric convention to be r}""^ = diag(l, ±1). 
A choice of gamma matrices vahd for both signatures is 

1 /I \ 2 fO -k 



^ \0 -l) ' ^ \k 

where k = i for Euclidean and k = 1 for Lorentzian signature. 

Using the gamma matrices and ([T|), we can rewrite Dirac's equations ^ as 

Adiip + Bd2->p + Cip- X-tp = 0, 

where 

-'l-U' '-i-l-'sy 



and 



Ai = ie\, A2 = —ike\, Bi = ief, B2 = —ike^, 



Let (x, y) be a coordinate system on the two-dimensional manifold and 



We now make the important assumption that the spinor is multiplicatively separable 

ipi = ai{x)bi{y). 

The Dirac equation then reads: 

Aiaibi + ^202^2 + Biaihi + -6202^*2 + Ciaibi - 6*20262 - Aoifti = 0, 

-A2aibi - Aia2b2 - -620161 - -810262 + 620161 - 610262 - A0262 = 0. (13) 

We now apply the general results of j20j to the Lorentzian case. 

Definition 1. The Dirac equation (jl3p and the operator D are separated in {x,y) if there exists 
nonzero functions Ri{x,y) such that (jlSp can be rewritten as 

RiaMEf + E'O = 0, R2aMEl + E^) = (14) 
for suitable indices r, s, t, u where Ef{x,aj,dj), E-{y,bj,bj). Moreover, the equations 

Et = ^^ = -Ey (15) 
define the separation constants //j. 
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The above definition refers to so-called "naive" separation of variables that is not always the 
most general. In order to find symmetry operators, we adopt for our analysis some assumptions. 

First of all, we assume a given coordinate system (x, y) and impose a separation of (jl3p 
according to the previous definition. This can be done in three different ways 

Type I: ai 7^ 02 and 61 7^ 62. 

Type II: ai = 02 and bi ^ 62 (or viceversa). 

Type III: oi = 002 = a and bi = db2 = b {c,d £ C). 

Following the procedure laid out in [20], we build eigenvalue-type operators = flip with 
eigenvalues ^(/Uj) making use of the terms Ef and in (jlSp only. We require that the opera- 
tors L are independent of A. Furthermore, A 7^ 0. 

Finally, we require the symmetry condition, that is [L, D]'0 = (LB) — DL)^/; = for all V- 
A operator L which satisfies the condition, is called a symmetry operator since it maps solutions 
into solutions. 

The symmetry operators are directly generated by the separated equations and having them 
enables one to immediately write down the same separated equations. In addition, the separation 
constants are associated with eigenvalues of symmetry operators. 

The only relevant case is Type I separation, since it is the only one associated with non- 
trivial second-order operators. We shall now make use of the naive separation assumption and 
of definition ([H]). We would like to determine the indices r, s, t, u and Ef{x, aj, dj), E^{y, bj, bj). 
Thus, we focus on what we can factorize from ()13p . 

By inspection, we notice that the only allowed factorizations are: 

• Factorize 0162 in the first equation and 0261 in the second equation. 

• Factorize 0261 in the first equation and 0162 in the second equation. 

We consider the first factorization, which implies Ai = B2 = 0. We divide into two parts both 
the equations: one part which is a function of x and the other a function of y, in order to 
have Ef, Ef and E^, E^. 
Hence 




Ef 




-Ml, 



and 




Ey 



i?li?2^+C,i?2^ + Ai?2^ 

bi oi oi 



where and fi2 are the separation constants. 

The second-order operator L is therefore defined by the following equations 



A2Rid2 - C2Ria2 = /^lai 



A2Ridi - C2Riai = /X2a2- 



If we set Hifi2 = /^i the operator is given by 






A2Rld^A2 - C2A2RI - C2RIA2 
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+ 



/ Cl{x,y)Rl{y) 
- A2{x,y)Rl{y)d^C2{x,y) 







V 







Cl{x,y)Rl{y) 
A2{x,y)Rl{y)d,C2ix,y)J 



We notice that the functions A2R1 and C2R1 are functions only of x. 

Finally, we look for the conditions that have to be applied in order to have the operator 
commuting with the Dirac operator D: it results already that [L,B] = 0, so no other conditions 
are needed. It follows from a detailed analysis of this case (D5 separation scheme in [20]) that 
we obtain a Liouville metric with one ignorable coordinate. This case will be discussed in detail 
in the next section. Another case is also possible which corresponds to the D7 separation scheme 
in [20] . However, it may be shown that it is equivalent to previous one modulo the sign of the 
Lorentz metric. 

Separability for equations of Type II and of Type III gives rise to first-order operators. 



5 Liouville metric 

In Section [2] we concluded our analysis of symmetry operators with the condition (|lip on the 
second-order operator: 

= -IVbiRe'^"). 

Its analysis requires knowledge regarding which spin manifolds M admit nontrivial valence two 
Killing tensors. 

To further proceed with our analysis, we recall the following important result. 

A two-dimensional Riemannian or Lorentzian space admits a non-trivial valence two Killing 
tensor if and only if it is a Liouville surface in which case there exists a system of coordinates 
(n, v) with respect to which the metric g and Killing tensor K have the following forms 

g = {A{u) + B{v)) {du^ + rjdv'^), 

B(v) ^ ^ A(u) ^ ^ 

A[u) + B(v) A{u) + B(v) 

where A and B are arbitrary smooth functions. Furthermore, the frame component of K are 
given by 

[X"''] =dmg{B{v),-TiA{u)). 

The spin frame corresponding to D5 separation discussed in the previous section is given by 





^A{u)+B{v) 



^A{u)+B{v) 





where ^ = and B{v) = l/Ri{v)'^. 

It follows from the previous section that the only Type I separation, other than the equivalent 
one discussed at the end of the previous section, is associated with the nonsingular Dirac operator 
and associated symmetry operator of the form 



Ri{y) 
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The operator K written above agrees with the second-order operator of Section [3] computed for 
the Liouville metric under the assumption a = = 0. We observe that it is in fact the square 
of the first-order operator 



corresponding to the Kilhng vector associated to the ignorable coordinate x of the Liouville 
metric where, as we said before, k = i for Euclidean and k = 1 for Lorentzian signature. 

The corresponding coordinates separate the geodesic Hamilton-Jacobi equation. If the mani- 
fold is the Euclidean plane, the coordinates, up to a rescaling, coincide with polar or Cartesian 
coordinates. In the Minkowski plane the coordinates correspond to pseudo-Euclidean or pseudo- 
polar coordinates. 

6 Separation in complex variables 

On real pseudo-Riemannian manifolds the Hamilton-Jacobi equation can be separated not only 
in standard separable coordinates but also in complex variables [8]. As in classical separation of 
variables theory, complex separable variables are determined by eigenvalues and eigenvectors of 
second-order Killing tensors. If the manifold is pseudo-Riemannian, pairs of complex-conjugate 
eigenvectors and eigenvalues of second-order real Killing tensors can exist in some part of the 
manifold, together with real ones. Where complex eigenvectors appear, it is impossible to 
determine real orthogonal separable coordinates, and the introduction of complex variables is 
necessary. The complex variables behave in all respects as complex coordinates, but they are 
not independent because of the conjugation relation. In the following their lack of independence 
will be irrelevant. Let us consider the 2-dimensional Minkowski manifold with pseudo-Cartesian 
coordinates {x,y). The geodesic Hamiltonian is given by 



The space of valence two Killing tensors is 6-dimensional and there are ten different types of 
separable webs real in some part of the space [7]. Another separable web, everywhere complex, 
is defined by [8] 



and is determined by the eigenvectors of the Killing tensor whose non-null components in (x, y) 
are 



and the associated polynomial first integral is 

L =PxPy 

By defining the canonical momenta as 
we have 





z = X + iy 



z = X — ly, 
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and a real complete separated integral of the Hamilton-Jacobi equation can be determined [8] . 
Because (z, z) are both ignorable variables, they should also separate the Dirac equation (indeed, 
in Minkowski space they are the only complex separable web with at least one ignorable variable) . 
With respect to (x, y) the Dirac operator can be written as 



.2 9y 



where (ea) are the components of the spin frame. Because the (ca) can be assumed to be 
constant, and since 

dz = ^{dx - idy), 
we have 



dz = ^{dx + id. 



yi-' 



e\ + iel 



'1 ^^1, 

In order to write D in the form 

1 or^^+Ho 



dz + 



-1 



corresponding to Type I separation in the Minkowski space, the components of (ca) in {x,y) 
must be 



2' 



1 

'2' 



el 



1 

2' 



el 



Therefore, by using (ca) and (eJJ) = (ea) ^ for raising and lowering indices. 



(ca) = lidx + idy, dx - idy). 

It is remarkable that the spin frame base allowing the separation of variables in (z, z) is essentially 
coincident with {dz,dz)- 

Both z and z are ignorable variables, therefore, both 



1 
1 
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and 



ij^' 



can be used as differential operators associated with separation of variables. The integration of 
the ipi = aihi is in all respects the same as for the separation in real coordinates. 

7 Appendix 

System ([7]) is obtained by using equations that characterize second- and third-order covariant 
derivatives. Such equations make use of Ricci's identities ([3]) and dH. 
For second-order covariant derivatives the following equation hold 

VaVcV' = |[Va, VelV^ + V aclp = ac'>P + f EacT^- 

This last equation can be rewritten as 



(16) 
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Similar equations for third-order covariant derivative require more calculations and we will show 
only the main passages 

Now by using (jl6p and expanding: 

Now the first part of the right hand side is just Vabc- The second part can be rewritten to obtain 

+ febc-fVa + ^{-riacVb - Vbc^a + 2r]abVc)- (17) 

Equations (jl7p and (|16p together with Ricci's identities are all the tools needed to expand 
equation ^ in system d?]). 
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